Dynamical fluctuations in an exactly solvable model of spin glasses. 
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In this work we calculate the dynamical fluctuations at 0(1/N) in the 
low temperature phase of the p = 2 spherical spin glass model. We study 
the large-times asymptotic regimes and we find, in a short time-differences 
regime, a fluctuation dissipation relation for the four-point correlation func- 
tions. This relation can be extended to the out of equilibrium regimes intro- 
ducing a function X t which, for large time t, as t~ x ^ 2 as in the case of the 
two-point functions. 



PACS Numbers 05-7510N 



I. INTRODUCTION 

The mean field Langevin dynamics for spin glasses has been quite extensively studied 
in recent years |!|, jj), j|, Q, (s) . The main result has been that for low enough 
temperatures there is an off-equilibrium regime where the dynamics of the system depends on 
its whole history up to the beginning of its observation and often this feature is accompanied 
by a loss of the validity of the fluctuation dissipation theorem (FDT) . Though it is in general 
very difficult to calculate the explicit time dependence of the correlation functions, except 
for some simple models, [||, their behaviour in all the different asymptotic regimes has been 
well understood || j| . 

The main question we are left with is how to extend the mean field picture for finite 
dimensional systems. From the analytical point of view this amounts to take into account 
the corrections to the mean field limit. 

In all generality one can determine closed equations for the two-point correlation func- 
tions as saddle points of an appropriate functional of a two time dependent field Q"^(ti, £2), 
0| . In this formalism, to consider the dynamical fluctuations around this mean field limit, 
one has to solve the equations of motion for the propagators, which are some four-times 
correlation functions and are related to the dynamical spin glass susceptibility. 

At the dynamical critical temperature, where dynamical scaling is supposed to hold 
and the off-equilibrium features are not relevant, these equations have been solved for some 
models M. 
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Unfortunately it is not clear how to approach the equations of motion below T c unless 
one knows which are the large-time asymptotic regimes for these four-time functions at low 
temperatures. 

Furthermore one would like to understand if, and in which time regimes, a kind of 
fluctuation-dissipation relation can be written for these four-time dependent quantities. 

These questions can by answered in the case of a simple spin-glass model which can be 
solved explicitly. It describes soft spins interacting through quenched random long range 
couplings and forced to satisfy a global spherical constraint. Statically, this so called spherical 
p = 2 spin glass does not present the peculiar features typical of other spin glasses and 
reveals itself to be a sort of disguised ferromagnet. Below a critical temperature T c the 
system freezes in one of two possible states spontaneously breaking the symmetry under the 
parity operator that connects them. 

Nevertheless, for temperatures below T c and if the dynamics of the system begins in 
a random initial configuration, the model presents a non trivial out of equilibrium dynam- 
ical regime exhibiting the so called aging phenomenon and a violation of the fluctuation- 
dissipation theorem. 

The simplicity of this model allowed for an explicit calculation of the two-point corre- 
lation and response functions In this work we shall calculate the explicit form of the 
four-point correlation functions. Hopefully, this simple case will provide us some hints on 
how to approach more complicated models which cannot be solved without the use of the 
functional methods and therefore without an ansatz on the asymptotic behaviour of the 
solution. 



II. THE STATICS 



The spin glass spherical model is described by the Hamiltonian 

H = --^JijSiSj -^hiSi, (1) 

where hi (t) is an external magnetic field, and the spin variables Si are forced to satisfy the 
spherical constraint 

N 

E S * = N - ( 2 ) 

i 

The couplings are symmetric quenched random variables extracted from a Gaussian 
probability distribution with Jy = and Jfj = 1/N. The mean field solution of the statics is 
rather simple for this model p0| . Note that, using conventional notations, we have indicated 
by (.) the averages over the Boltzmann distribution, and by t the averages over the different 
realizations of the quenched disorder Jy . 

The partition function is 
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where we have set hi = for all i and z is a Lagrange multiplier introduced to enforce 
the constraint. For N — > oo the distribution of the eigenvalues of the matrix follows the 
Wigner semi-circle law H]: 

P(A) = ^V4-A2 , Ae[-2,2]. (4) 

In this limit one can formally write the saddle point equations on Z which reproduce a softer 
version of the constraint 

, N f2 1 

In equation (Q) s a = Y2i=i < l ) ? s i * s ^ ne projection of the spin variables on the a-th 
eigenvector of the matrix of the couplings . 

It can be seen that below a critical temperature T c = 1 the second of the equalities in 
^ does not hold because a spontaneous magnetization arises along the eigenvector with 
eigenvalue 2 (or —2). For T < T c the value of the Lagrangian multiplier remains fixed at 
the branch point z = (3 and a magnetization (S2) oc N 1 / 2 appears. 

The spin glass susceptibility is defined as follows: 

ij a ^ ' 

In the large- iV limit XSG diverges as 1/(T — T c ) for T — > T c + and it remains infinite in 
the whole frozen phase. This model is critical for all temperatures below T c . 

The computation of the average free energy can also be done using the replica method, 
H , where a replica symmetric ansatz solve the model exactly. One can reproduce the results 
for Xsg using the following identity: 

XSG = lim [(Sqlp) - 2{5q a p8q ai ) + (5q a p5q 7 s)] , (7) 

where the 5q a p are the fluctuations around the saddle point of the replicated partition 
function. It can be seen that each of the terms of equation (m) diverges below T c . 



III. THE MEAN FIELD DYNAMICS 

We shall mainly deal with the dynamical behaviour of the model. A Langevin dynamics 
for the Hamiltonian ([I]) joined with the spherical constraint gives the following equation of 
motion 

^ = E - w*) - *(*)*(*) + m*) + &(*) (8) 

3 

where z(t) is a time-dependent Lagrange multiplier, and is a Gaussian noise with zero 
mean and variance (^i(t)^j(t r )) = 2TSijS(t — t'). In this model the solution of (||) for the 
component s a can be explicitly written 
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f.\ «(*-*o)-f* z(t)cIt f a {t-t")-V z{t')<It' . „. ,j, u ,,„ ,„x 

where again a labels the eigenvectors of Jy, and the initial time is t = 0. We introduce the 
correlation and response function which read: 



1 N 

C(t 1 ,t 2 ) = -^( Si (t 1 ) Si (t 2 )), (10) 

i 

n(+ +\ 1 djsjjh)) 1 A djsijh)) 1 ^ 1 

Gfo, fc) = ^ r = ^ "flFTTX =«L ~f <*(*i)6(*a)>, (11) 



where the last equality is valid because the noise has a Gaussian distribution probability. In 
the N — > oo limit and in absence of the external field the above functions have the form S 



C(h,t 2 ) 



1 



rt2 

2T / dt'T(t 
Jo 



(h + h) J y ' {tx + h-2t>) 



(12) 



where I\ [x] is the modified Bessel function and the function 

2 f z(t)cLt 

r(t) = e J *o v j (14) 



is fixed by implementing the spherical condition C (t, f ) = 1 in equation (|l£ 

The simplicity of this model allows one to obtain the mean field solutions ([l2]) ( |l3| ) 

directly, simply averaging over the distribution of the eigenvalues of the matrix Jij . 

For other models of spin glasses we do not manage to obtain the explicit form of C(t\ , < 2 ) 

and G(t%, i 2 ). 

A quite general procedure, introduced in jy, allows to obtain closed equations for 
C{t\,t2) and G(ti,f 2 ) as the saddle point solutions of a dynamical generating functional. 
In particular, in the case the off-equilibrium dynamics of the p-spin model, with p > 2, 
and the SK model, we are not able to solve these equations and therefore to determine the 
explicit time dependence of the correlation and response functions but we can only predict 
the structure of their asymptotic behaviour. This model that provides us with an explicit 
solution, allows us to control the assumptions and the ansatz we used to obtain analytical 
results for the dynamics of more complicated systems. 

For the correlation functions C(£i,£ 2 ) and G(ii,i 2 ) it has been shown that there are 
basically two asymptotic time scales and that they can be distinguished by the variable 

\ = t 2 /h§. 

• For A ~ 1 the system is in an equilibrium regime, in which the functions depend only 
on the difference of the two arguments (time translational invariance) and the FDT relation 

G{t l -t 2 ) = l 9C{ \- h) e{t l -t 2 ), (is) 

holds at all temperatures. 
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• For A~0(1)<1 and for T < T c , the system is in the so called aging regime: the 
correlation functions depend on both time variables (in this case through A) and not only 
on time differences. Moreover the FDT relation is not valid although it can be generalized 
introducing a function X tl (C) || 0] by the relation 

G(hM = ^«hM m - h , (16) 

For large times one has that X tl (C) — > X(C). The function X(C) generally characterizes 
the type of aging dynamics of the model. In this model the slowness of the dynamics is 
due to the flatness of the energy landscape and is therefore qualitatively similar to ordinary 
domain coarsening. This kind of dynamics is in general associated with a X(C) = which is 
indeed the case for this model in which it has been found that, in the aging regime, Xt x (C) 
scales as t^ 1 ^ 2 for large t%, ||. We shall see that this will be the case also for the four-times 
functions. 



IV. DYNAMICAL FLUCTUATIONS 

The problem of the dynamical fluctuations around the mean field solution for spin glass 
models, for T <T C has not been faced yet. We shall therefore study the four-point functions 
in this model in which it is possible to determine the explicit temporal behaviour of these 
functions and we will describe the different large-time asymptotic regimes. 

Let us now introduce the four-time correlation functions 

lihMMM) = ^ J2 (WOT^W ~ (*(*i)si(t 2 )> (0fe)0(i 4 )>) (17) 
KtiMMM) = ^ J2 - (si(*i)si(t 2 )> (sjfo&fa))) (18) 

((siW-SifamiisjW-sjiU))*)). (19) 

The functions 7, A, u) are related to the fluctuations around the saddle point of 
Qf^(ti,t2), [0 that gives, as mean field solutions, the two-point correlations and re- 
sponse functions. Note that 7(^1^2,^3,^4) is related to the dynamical XSG'- 

= 1 v djsjjh)) gggy 

N% dh^ts) dh 3 (t A ) ' 1 ! 

while uj{t\, t2, £3, £4) is the dynamical four-points correlation function defined in such way to 
stay finite in the asymptotic limit of t\ — > 00 with t\ ~ t% ~ £3 ~ t&, i.e. on the equilibrium 
time scale. 

In this regime one expects to verify some FDT-like relations 
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-^j-L0(tl,t2, t 3 , ti) — 2T^2A(tl,t2, h, ti) — \(tl,tl,t 3 , ti) — X(t 2 ,t2, t 3 , t±) 

—2\(ti,t 2 ,t i , i 4 ) + X(ti, tift^, t/j) + \(t 2 ,t 2 ,t i , i 4 )^ , 

d 

— -X(ti,t2,t3,t 4 ) — Tj(ti,t2,t 3 ,t4). 
Ot 3 



(21) 
(22) 



Let us remark that in equation (|17|) (|18|) (19) we have by hand subtracted the part of 
the term that is of order O(l), and the quantities defined are all of order 0(1/N). 

In the functional formalism the above functions are the propagators of the fields 
Q a ^{t\,t2) introduced in [jl| evaluated at zero momentum. This simple model will pro- 
vide us with the first explicit calculation of the dynamical propagators for a spin glass in 
the low temperature phase. 

Let us now set the time order t± < t 3 < £2 < t%. Using the fact that the model is 
quadratic we can write 



7(*1>*2,*3,*4) = 7(£l,i2,£3,i4) +7(*l!*2,*4,*3) 

A(ti, t 3 , £4) = X(ti, t 2 , t 3 , £4) + \(t 2 , ti, £3, £4) 

W(ti,t 2 ,*3,*4) = 4w(il,£2,*3,*4) + -Hj(ti,t 2 ,t 4 ,t 3 ) - 4w(ii,ti,t3,t 4 ) 

4L)(t 2 ,t 2 ,t 3 ,t 4 ) - 4u)(ti,t 2 ,t 3 ,t 3 ) - 4w(h,t2,ti,U) + 
2uj(t 1 ,t 1 ,t 3 ,t 3 ) + 2Q(t2,t2,t 3 ,t 3 ) + 2w(t 1 ,t 1 ,t A ,U)+2u(t2,t2,U,U) 

where we have defined 



lihMMM) = {s i {t 1 )^{t 3 )){s i {t2)^{U)) 



X(ti,t2,t 3 ,U) = (Si(ti)Sj(t 3 )){Si(t 2 )£j(t4,)) 



Uj{ti,t 2 ,t 3 ,t i ) = (Si(ti)Sj(t 3 )}{Si(t2)Sj(t 4 )}. 



One gets 

7(*1)*2,*3,*4) 

X(ti, t 2 , t 3 ,t 4 ) 

U)(tl, t2, t 3 ,t 4 ) 



—e{ti-t 3 )0{t2-u)\ 



n ^r(< 1 )r(t 2 )r(i 3 



lT(t 3 )T{t 4 ) I^jh + t 2 - t 3 - t 4 )} 

T(h)T(t2) (h + 1 2 - h - u) 

h[2{t 1 +t 2 +t 3 ~U)} 



(h +t 2 + t 3 - U) 



-2T 



1 



t3 dt , r(t , ) h[2(t 1 +t 2 +t 3 -t 4 -2t') 
[ 1 (ti+t 2 +t 3 -t 4 -2t') 
1 



4T" 



N V / 1)WW 

Jo {t 1 +t 2 +t 3 +U-2t') 

2T [\'T(t') h[2{h+h + h + t '- 2t ' )] + 
Jo {t 1 +t 2 +t 3 + t 4 -2t') 

h [2(t! +t 2 +t 3 + t A - 2t' - 2t")] 



dt' dt"T(t')T(t") 



{ti +t 2 + t 3 +U- 2V - 2t") 



(23) 
(24) 



(25) 

(26) 
(27) 
(28) 

(29) 



(30) 



(31) 
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We shall study the fluctuations around the mean field solutions and in particular we 
shall be interested in understanding if equation s (pl|) a n d are verified. 
In the large-time limit form the functions (|29|) (|30|) (|3l|) read 



j(tx,t 2 ,ts,t 4l ) = 
A(ti, t 2 , ti) = 



ht 2 



-2[*i+ia-ts-t4] h[2(h +t 2 -t 3 - U)] 

(h + t 2 -t 3 - U) 



(32) 



4(*1*2* 3 ) 



U(h + t 2 + t 3 -U) 2 _ 

2{t 1 +t 2 -t 3 -t 4 ) 

1 - T I dt'- 



h[t'} 



f(l 



2{t 1 +t 2 +t 3 -t i ) ' 



(33) 



u){ti,t 2 ,t 3 ,U)= q EA V4:ir 



dt" 



T(t")e~ 4t " 



(1 



21'' 



tl+*2+i3+*4 



)3/2 



(tl + *2 + *3 + i 4 ) 2 

■2(t 1 +t 2 +t3-*4-2t") 



[2(l-q 2 EA )-2T + 2Tx 



1-Qea-T 



dt'- 



[t'(l 



2(t 1 +t 3 +t 2 +t 4 -2t") 



ttt) 3 / 2 



(34) 



3 /2 

Note that the function uj is divergent as in all the time regimes. This is due to the fact 
that w(ii, t2 ; ^3) ^4) diverges in the static limit. 

However, in the definition (19) of the u>(ti,t 2 ,t3,t4) we have subtract this overall be- 
haviour on the equilibrium time scale, and in this regime we can consider only the dynamics 
on the time differences. 

The structure of the possible time regimes is in principle similar to that for the two 
point correlation functions: the times are either far or close to each other. 

In a generic four-time function we can have some times which are close and some which 
are far away from each other and this complicates the separation in time sectors of the time 
space because the function may be in local equilibrium with respect to some times and aging 
with respect to others. We will now study the functions in all the possible times regimes 
obtainable by four ordered times. We will calculate their asymptotic behaviour for each 
regime. The result of our analysis is that relations (|2l]) and ( |22| ) hold if and only if all the 
four times are close to each other. That means if there are no aging time scales. Otherwise 
the dependence upon the aging times (i.e. they are far from each other) will dominate the 
functions. 

The results that we obtain are summarized in the table below where with we 
indicated in a compact form the r.h.s. of the equation (21). On the first column we indicated 
the time regime we were considering, using curly brackets to group times whose difference 
remains finite in the limit of t\ — > 00. 
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1 imeKegime 


U! 


U 


2^ A 


A 


d\ 
*3 


7 




,-5/2 
*1 




0(1) 


0(1) 


r 7 


0(1) 


(flt 2 )(*3*4) 


,-5/2 
*1 


,-5/2 
*1 


, — 3 

h 


0(1) 


,-1 
*i 


,-3/2 
h 


11^2 (.C3t4 J 


,-1/2 
*1 


,-1/2 
*1 


,—1 

h 


0(1) 


,—i 

k 


,-3/2 
k 


, [ii Si 


,3/2 
*1 


,1/2 
*1 


0(1) 


0(1) 


h 


,-3/2 
*i 




,-1/2 


,-3/2 
*1 


+-2 


0(1) 


h 


,-3/2 
*i 


{t 1 £2^3)^4 


,-1/2 
c l 


,-3/2 
t l 




0(1) 




1-3/2 
n 


^lfe^M 


-1/2 
l l 


-1/2 
c l 




0(1) 




1-3/2 
l i 


t^i^ti 


,3/2 
l l 


1/2 
E l 


0(1) 


0(1) 


t-, 1 


1-3/2 



To see if equations (^l|) and (|2^) hold one has to compare the third column with the 
fourth and the sixth column with the seventh. One sees that the FDT relations for the 
fluctuations hold only in the first regime, where all the times are at a finite distance from 
each other. In the other regimes they still can be generalized introducing an X tl (C) oc t 1 1 
as in the case of the two-times functions. Let us also remark that u>(ti, £2, £3, £4) diverges 
every time t\ — ti and £3 — £4 are order 0(t\). In this cases, the function u)(ti, £2, £3, ti), 
(19), has not a well defined limit. Moreover, in the FDT regime the functions 7(^1, £2, £3, £4), 
A(ii, £2, £3, £4) and u)(ti, £2, £3, £4) are asymptotically dependent on the times only trough 
t = ti + ti — £3 — £4, while in the aging regimes an explicit dependence from the times that 
are at distance O(ii) remains. Finally we would like to point out that r is a characteristic 
time scale that determines whether we are in an FDT regime (r ~ O(l)) or not (r ~ 0(ti)), 
assuming the same role as t — t' in the case of the two times functions. 



V. CONCLUSIONS 



In this work we calculated explicitly the four-times functions for the spherical 2-spin 
glass model in the low temperature phase. We calculated the asymptotic behaviour of these 
functions in the possible time regimes that can be selected with four ordered times. We 
found out that there is a time regime in which one can relate the functions through a relation 
similar to the fluctuation-dissipation theorem for the equilibrium dynamics of the two-times 
functions. In this work we provide the first explicit expression for the dynamical four-point 
functions for a spin glass model in the cold phase. It is our intention to use the results 
obtained here to analyze less trivial spin glass models describing short range interactions of 
p > 2 spins. In these models there is lot to be understood about the dynamical propagators 
and about their static limit . Unfortunately the equations for the propagators are integro- 
differential equations that cannot be solved without the use of an ansatz on their asymptotic 
form. We think that the knowledge of the dynamical structure of the propagators for the 
spherical p — 2 spin glass model can help us in view of these further developments. 
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